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q ■ Abstract 

We give a classical protocol to exactly simulate quantum correlations implied by 
a spin-s singlet state for the infinite sequence of spins satisfying (2s + 1) = 2 n , in 
' the worst-case scenario, where n is a positive integer. The class of measurements 

\f} ■ we consider here are only those corresponding to spin observables. The required 

amount of communication is found to be logfcd where d = 2s + 1 is the dimension 
q | of the spin-s Hilbert space. 

\Q '. 
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3 1 Introduction 

cr 

. !_( ! It is well known that quantum correlations implied by an entangled quantum state of a 

bipartite quantum system cannot be produced classically, i.e., using only the local and 
realistic properties of the subsystems, without any communication between the two sub- 
systems [1]. By quantum correlations we mean the statistical correlations between the 
outputs of measurements independently carried out on each of the two entangled parts. 
Naturally, the question arises as to the minimum amount of classical communication 
(number of cbits) necessary to simulate the quantum correlations of an entangled bipar- 
tite system. This amount of communication quantifies the nonlocality of the entangled 
bipartite quantum system. It also helps us gauge [2] the amount of information hidden in 
the entangled quantum system itself in some sense, the amount of information that must 
be space-like transmitted, in a local hidden variable model, in order for nature to account 
for the excess quantum correlations. 
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In this scenario, Alice and Bob try and output a and (3 respectively, through a classical 
protocol, with the same probability distribution as if they shared the bipartite entangled 
system and each measured his or her part of the system according to a given random 
Von Neumann measurement. As we have mentioned above, such a protocol must involve 
communication between Alice and Bob, who generally share finite or infinite number of 
random variables. The amount of communication is quantified [3] either as the average 
number of cbits C(P) over the directions along which the spin components are mea- 
sured (average or expected communication) or the worst case communication, which is 
the maximum amount of communication C W (P) exchanged between Alice and Bob in any 
particular execution of the protocol. The third method is asymptotic communication i.e., 
the limit lim n ^ OQ C(P n ) where P n is the probability distribution obtained when n runs 
of the protocol carried out in parallel i.e., when the parties receive n inputs and produce 
n outputs in one go. Note that, naively, Alice can just tell Bob the direction of her mea- 
surement to get an exact classical simulation, but this corresponds to an infinite amount 
of communication, the question whether a simulation can be done with finite amount 
of communication was raised independently by Maudlin [4], Brassard, Cleve and Tapp 
[5] and Steiner [6]. Brassard, Cleve and Tapp used the worst case communication cost 
while Steiner used the average. Steiner's model is weaker as the amount of communica- 
tion in the worst case can be unbounded although such cases occur with zero probability. 
Brassard, Cleve and Tapp gave a protocol to simulate entanglement in a singlet state 
(i.e., the EPR pair) using eight cbits of communication. Csirik [7] has improved it where 
one requires six bits of communication. Toner and Bacon [8] gave a protocol to simulate 
two-qubit singlet state entanglement using only one cbit of communication. Interestingly, 
quantum correlations that cannot be classically simulated without communication also 
occur in a scenario where incompatible observables are successively measured on class of 
input (single particle) spin-s states which can be simulated with a classical protocol with 
communication between successive measurements [9]. 

Classical simulation of quantum correlations is accomplished for spin-1/2 singlet state, 
requiring the optimal amount, namely, 1 cbit of classical communication in the worst-case 
scenario, using arbitrary projective measurement on each site [8]. It is important to know 
how does the amount of this classical communication change with the change in the value 
of the spin s, in order to quantify the advantage offered by quantum communication 
over the classical one. Further, this communication cost quantifies, in terms of classical 
resources, the variation of the nonlocal character of quantum correlations with spin values. 
In this paper we give a classical protocol to simulate the measurement correlation in a 
singlet state of two spin-s systems, considering only measurement of spin observables 
(i.e., measurement of observables of the form a. A where a is any unit vector in R and 
A = (A x , A y , A z ) with each A$ being a (2s + 1) x (2s + 1) traceless Hermitian matrix and 
the all three together form the £77(2) algebra) where 2s + 1 = 2 n , n being any positive 
integer. As s can only take integral or half-integral values, the allowed values of s form 
the infinite sequence s n= i = 1/2, s n=2 = 3/2, s n=3 = 7/2, etc. For s n = 2 n ~ 1 — 1/2, our 



2 



protocol requires n cbits of communication and 2n number of independent and uniformly 
distributed shared random variables. 

We describe the general classical simulation of singlet state of two spin-s systems in 
section II, followed by the specification of the measurement scenario, which we will be 
considering in the present paper. We briefly describe the protocol of Toner and Bacon 
[8] in section III. We present our results for spin-s singlet state in section IV. Finally, our 
conclusions are summarized in section V. 

2 Correlation of two spin-S* singlet state 

The singlet state \i/j~)ab of two spin-s particles A and B is the eigenstate corresponding 
to the eigenvalue of the total spin observable of these two spin systems, namely the 
state 

\1>7)ab = — L= £ (-l) s - m \m) A ® I - m) B , (1) 

where | — s), | — s + 1), . . ., \s — 1), \s) are eigenstates of the spin observable of each of 
the individual spin-s system. Thus ab is a maximally entangled state of the bipartite 
system A + B, described by the Hilbert space C 2s+1 <8> C 2s+1 . In the case of classical 
simulation of the quantum correlation (i/j~, 2 \a.a <g> b.a\ip^ 2 ) of the two-qubit singlet state 

\ipy 2 ), Alice considers measurement of traceless observable a. a and Bob considers that 

of the traceless observable b.a . These are spin observables. Analogous to the Pauli 
matrices, one can consider (2s+l) 2 — 1 number of trace-less but trace-orthogonal Hermitian 
(2s + l) x (2s+l) matrices Ai, A 2 , . . ., A( 2s +i)2_i (i.e., TrA; = for all i but Tr(AjAj) = if 
i 7^ j; see, for example, [10] ) such that a general projective measurement on the individual 
spin-s system corresponds to the measurement of an observable of the form c.A, where c 
is a unit vector in 

R (2s+if-i and £ ig the + ^ 2 _ !)_ tuple (Ax, A 2 , . . . , A (2s+1)2 _ 1 ) of 

the above-mentioned A matrices. In general, the quantum correlation (ip~\c.A ® d.A\ip~) 
will be a bilinear function in the components on c and d. In the special case when 
(Aj ® Aj)\ip~) = aSij\^~) for all i,j — 1, 2, . . . , (2s + l) 2 — 1, the quantum correlation will 
be of the form ac.d . 

Classical simulation of this general quantum correlation seems to be quite hard one 
possible reason being the absence of Bloch sphere structure for higher spin systems. 
Rather we will consider only measurement of spin observables, namely the observables 
of the form a. J on each individual spin-s system, where a is an arbitrary unit vector in 
R 3 and J = (J x , J y , J z ). For the (2s + 1) x (2s + 1) matrix representations of the spin 
observables J x , J y , and J z , please see page 191 - 192 of ref. [11]. 

J matrices satisfy the SU(2) algebra, namely [J x , J y ] = iJ z , [J y , J z ] = iJ x , [J z , J x ] = 
iJ y . The eigenvalues of a. J are — s, — s + 1, . . ., s — 1, s for all a G -R 3 . The quantum 
correlations (ip~\d.J ® b.3\ip~) (which we will denote here as (a/?), where a runs through 
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all the eigenvalues of a. J and (3 runs through all the eigenvalues of b.J) is given by 

<^-|a.J <g> 6.J|VJ> = = + l)a.b , (2) 

where a and b are the unit vectors specifying the directions along which the spin compo- 
nents are measured by Alice and Bob respectively [12]. Note that, by virtue of being a 
singlet state ,(a) = = ((3) irrespective of directions a and b. From now onward, we will 
consider only those spin-s systems for which 2s+l = 2 n , n being any positive integer. Thus 
we see that the allowed spin systems will form the sub-class {2 n_1 — 1/2 : n — 1, 2, . . .} of 
half-integral spins. 



3 Classical simulation of two spin- 1/2 singlet state 

As the working principles of our protocol are of similar in nature with those of Toner and 
Bacon [8], before describing our protocol, we would like to briefly describe the protocol 
of Toner and Bacon to simulate the measurement correlations on \ip^j 2 )- I n this scenario, 

Alice and Bob's job is to simulate the quantum correlation {ip~, 2 \d.^a®b.^a\^, 2 ) = — jd.b, 
together with the conditions that (a) = = ((3). To start with, Alice and Bob share two 
independent random variables A and £i, each of which has uniform distribution on the 
surface of the Bloch sphere S 2 in R 3 . Given the measurement direction a, Alice calculates 
— 7;sgn(a.\), which she takes as her measurement output a. Note that sgn(x) = 1 for 
all x > and sgn(x) = — 1 for all x < 0. As A is uniformly distributed on S2, for 
each given a, — \sgn(a.\) will take its values | and — | with equal probabilities, i.e., 
Prob(o; = 1/2) = Prob(o; = —1/2) = 1/2 (and hence, (a) = 0). Alice then sends 
the one bit information c = sgn(a.\)sgn{a.fi) to Bob. Note that instead of sending 
sgn(a.\), by sending c, Alice does not allow Bob to extract any information about her 
output a. This is so because Prob(o; = l/2|c = 1) = Prob(a = — l/2|c = 1) and 
Prob(a = l/2|c = —1) = Prob(a = — 1/2 |c = —1). After receiving c, and using his 
measurement direction b, Bob now calculates his output (5 = ^sgn[b.(X + cfi)}. Now 

(P) = TuTTo I I sgn[b.(X + sgn(a.\)sgn(a.fi)fi)}d\dfi . (3) 
2 W J\es 2 Jfies 2 

Given any (1 £ S2, for each choice of A £ S2, the two values of the integrand corresponding 
to A and —A are negative of each other. As the distribution of A on S2 is taken to be 
uniform, the above-mentioned observation immediately shows that {(3) =0. As (3 £ 
{1/2,-1/2}, therefore Prob(/3 = 1/2) = Prob(/3 = -1/2) = 1/2. In order to compute 
(a/3), one should observe that Bob's output can also be written as (3 = ^X)d=±i[(l + 
cd)/2]sgn[b.(X + d p)]. The following two among the four integrals (which appears in 
(a(3)) ' 

- /, x 2 / / sgn(d.X)sgn[b.(X± fi)}dXdjl 
8(4vr) 2 J Xe s 2 J^s 2 
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cancels each other by incorporating the inversion jl — > —fx. And the rest two integrals 

± Q ,} x 2 / / sgn(a.X)sgn[b.(X± jl)}dXdfi 
8(4tt) 2 Jx eS2 Jixes 2 

are same and they are equal to the integral 

1 



sgn(a.\)sgn[b.(ji — X)]dXdjl . 

AeS 2 J fi£S 2 



8(4tt) 2 
And hence we have 

(ipy 2 \d.a <g) b.a\ipy 2 ) = (a/3) = 

/ / sgn(d.X)sgn[b.(jl— X)]dXdfi = — - / sgn(d.X)b.XdX = — d.b 

J\eS 2 Jfies 2 8(4tt) J~ X£ s 2 ' 4 

(4) 



8(4tt) 



4 Classical simulation of two spin-<S singlet state us- 
ing spin measurements 

Let us now come to our protocol. In the simulation of the measurement of the observable 
a. J (where a G R 3 is the supplied direction of measurement), Alice will have to reproduce 
the 2™ number of outcomes a = T~ x - 1/2, 2 n ~ l - 3/2, . . . , -2"" 1 + 1/2 with equal 
probability. If we consider the series —\ Ylk=i f(k)2 n ~ k , where, for each k, f(k) can 
be either 1 or —1, it turns out that the series can only take the above-mentioned 2 n 
different values of a. The probability distribution of these different values of the series 
will depend on that of the n-tuple {/(l), /(2), . . . , f(n)}. In order to make this probability 
distribution an uniform one (which is essential here for the simulation purpose), we choose 
here f(k) = sgn(d.Xk) for each k, where a is the measurement direction for Alice while Ai, 
A 2 , • • ., A„ are independent and uniformly distributed random variables on S2. We have 
seen in the above-mentioned Toner and Bacon protocol that if Alice and Bob share the 
two independent and uniformly distributed random variables A& G S2 and fik G S2, then 
the random variable r\ = sgn[b.(Xk + sgn(d.Xk)sgn(d.fik)fik)} is uniformly distributed over 
{1, —1}. Hence, as above, the quantity | Xlfc=i 2 n ~ k sgn[b.(Xk+ri : £ik)] will have 2 n different 
values = 2™" 1 - 1/2, 2 n ~ l - 3/2, ... , -2 n ~ 1 + 1/2 all with equal probabilities. But the 
interesting point to note is that in the calculation of the average (over the independent 
but uniformly distributed random variables Ai, A 2 , . . ., A n , p>i, £12, ■ ■ ■, An) of the product 
a/3, there will be no contribution from cross terms like sgn(d.Xk)sgn[b.(Xi +rifii)] if k ^ I. 

The protocol proceeds as follows: Alice outputs a = — \ J2k=i ^ n ~ k sgn(d.X k ). Alice 
sends n cbits Ci,C2, ■ ■ ■ ,c n to Bob where Ck = sgn{d.Xk)sgn{d.(xk) for k = 1,2, ... ,n, 
where Ai, A 2 , . . ., A n , fi\, £12, ■ ■ ■, fan are independent shared random variables between 
Alice and Bob, each being uniformly distributed on S 2 - Thus we see that, in terms 



5 



of shared randomness, A = (Ai, A2, . . . , A n , fa, fa, ■ ■ ■ , fa) is the shared random variable 
between Alice and Bob. After receiving these n cbits from Alice, Bob outputs (3 = 
\ YTk=\ 2 n ~ k sgn[b-(\ k + c k fa)}. It follows immediately from the discussion in the last 
paragraph that 

1 1 " f 

(a/3) = _____ ^2 2 2n ~ 2k J d\i... d\ k -id\ k+1 . . . dX n dfa . . . dfa-idfa+i . . . dfax 

/ sgn(d.X k )sgn[b.(X k + d k fa)]dX k dfa. (5) 

It follows from the discussion in section III regarding Toner and Bacon's work that (a (3) = 
~\ Z~2k=i 2 2n ~ 2k d-b ■ Summing the geometric series and using (2s + 1) = 2 n we finally get 

(af3) = --s(s + l)a.b . (6) 
3 

This protocol exactly simulates quantum mechanical probability distribution for par- 
ticular types of projective measurements, namely the spin measurement, on the spin s 
singlet state with 2s+ 1 = 2™ for positive integer n. The above protocol applies to infinite, 
although sparse, subset of the set of all spins (i.e., all integral and half integral values). 
The most important finding is that the amount of communication goes as Zog 2 (2s + 1) or 
as log 2 s for s ^> 1. Our protocol works equally for any two spin-s maximally entangled 
state as that can be locally unitarily connected to the singlet state. 



5 Conclusion 

Our result provides the amount of classical communication in the worst case scenario if 
we consider only measurement of spin observables on both sides of a two spin-s singlet 
state with the restriction that the dimension 2s + 1 of each subsystem must be a positive 
integral power of 2, and just n = log2(2s + 1) bits of communication from Alice to Bob is 
sufficient. We are unable to show whether our protocol is optimal (in the sense of using 
minimum number of classical communication). On the other hand, if we consider most 
general projective measurements on both the sides of a maximally entangled state of two 
qudits, with d = 2 n , it is known that (see [5]) Alice would require at least of the order of 
2 n bits of communication to be sent to Bob, in the worst case scenario when n is large 
enough. But for general d, log 2 d can be shown to be a lower bound on the average amount 
of classical communication that one would require to simulate the maximally entangled 
correlation of two qudits considering most general type of projective measurements [13]. 
It is also known that log 2 d bits of classical communication on average is sufficient to 
simulate the measurement correlation of a maximally entangled state of two qudits, when 
both Alice and Bob consider only measurement of traceless binary observables [14]. It 
thus seems that even if simulation of maximally entangled correlation in the most general 
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case of projective measurement is a hard problem, and one would require to send classical 
communication at least of the order of the dimension (for large dimensional case), there 
is still some room to search for efficient simulation protocols in lower dimensions. 
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